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Abstract 

Working in the framework of the ab —initio no-core shell model, we de¬ 
rive two-body effective interactions microscopically for specific harmonic- 
oscillator basis spaces from the realistic Argonne V8' nucleon-nucleon po¬ 
tential. However, our two-body effective interaction is only the dominant 
part of the A—body (A > 3) effective interaction, which also contains 
three-, four- and, in general, up to A—body components. We model 
these multi — body components as a simple two-parameter, basis-space- 
dependent interaction. We then investigate the coupling needed to pro¬ 
vide accurate descriptions of selected nuclear observables. We successfully 
model the multi-body interaction as a zero-range two-body operator times 
the number operator N sps , that measures the harmonic-oscillator quanta 
of the spectator nucleons. 


1 Introduction 

Despite the tremendous progress realized in the study of the bare nucleon- 
nucleon (NN) interaction as well as equally impressive advances in the theory 
of the effective NN interaction in the shell model, it has not been possible to 
systematically account for the properties of finite nuclei solely on the basis of 
two-body interactions. Studies of the roles played by many-body forces ( e.g . 3- 
body and 4-body) suggest that they provide small yet important contributions 
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[jl] U. However, more careful analysis of nuclear many-body systems is needed to 
fully identify these roles [T Adding to this complexity, the ab — initio No Core 
Shell Model ( NCSM ) [Q, |(]] generates effective many-body forces which have 

been shown to decrease in strength with increasing basis space size. Overall, 
the nature of the ’’real” and the ’’effective” many-body forces remains quite 
complicated and elusive. 

In this work, we present a phenomenological model of three-body and higher 
multi-body interactions as a simple, two-parameter and basis-space-dependent 
interaction. This multi-body interaction is added to our theoretically-derived 
two-body effective interaction, and its two parameters are adjusted in each basis 
space so as to reproduce results in agreement with experimental data. Thus, ours 
is a residual multi-body interaction obtained in a phenomenological way. Our 
chosen functional form is a delta function reminiscent of the recent attention of 
Effective Field Theory (EFT) to expansions of nuclear interactions beginning 
with zero-range terms. In EFT there is an appropriate cutoff in momentum 
transfers, corresponding to the poorly-known short-range part of the nuclear 
potential [JtJ. In a similar spirit, we model the poorly known many-body forces 
as harmonic-oscillator-cutoff-dependent zero-range interactions ||. 

More specifically, we model the residual multi-body interaction by a one- 
parameter two-body delta interaction and a one-parameter multi-body interac¬ 
tion. The multi-body interaction consists of a two-body delta times a number 
operator, N sps , that measures the excitations in oscillator quanta of the spec¬ 
tators. We expect and find that the contribution of the two-body part of the 
interaction weakens as the basis space increases and as our theoretically-derived 
effective two-body interaction approaches the exact two-body interaction. On 
the other hand, we find as expected that the multi-body part of our residual 
interaction remains relatively strong, which provides an indication of the role 
played by n —body forces (n > 2) in nuclei. 

2 The Shell-Model Hamiltonian and the Effec¬ 
tive Interaction 

The ab — initio approach in shell-model studies of the nuclear many-body prob¬ 
lem SII begins with the one- and two-body Hamiltonian for the A-nucleon 
system, i.e., 


i=l i<j 

with to the nucleon mass and V,j the NN interaction. We then add the center- 
of-mass harmonic-oscillator (HO) potential \m^L 2 R 2 , with 
R = (1/A) Yld=i O- The total center-of-mass harmonic-oscillator hamiltonian, 
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H^ m = P 2 m /2Am + iAmfi 2 P 2 , P c .m. = Ylt=\Pii later subtracted to leave 
a purely intrinsic Hamiltonian. At this stage we have a convenient mean field 
acting on each nucleon that helps our overall convergence when working in a HO 
set of basis states. We emphasize that we adopt proper measures to ensure that 
the intrinsic properties of the many-body system are not affected by the center- 
of-mass term and we discuss our specific procedure for this in more detail below. 
The modified Hamiltonian, thus, acquires a dependence on the HO frequency 
fl, and can then be written as 


hQ = E 


tr * | (~\2 

- —mil r, 

2m 2 
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E 

i<j 
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m£l 2 



( 2 ) 


Actual shell-model calculations can only be carried out in a finite basis space 
defined by a projection operator P, with the complementary space to the finite 
basis space (i.e. the excluded space) defined by the projection operator Q = 
1 — P. Furthermore, due to its strong short-range part, the realistic nuclear 
interaction in Eqs. (1) and (2) will yield pathological results unless we derive a 
basis-space dependent effective Hamiltonian: 


H. 


eff 
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(3) 

The second term in Eq. (3) is an effective interaction that is, in general, an A- 
body interaction. If the effective interaction is determined without any approxi¬ 
mations, the basis-space Hamiltonian provides an identical description of a sub¬ 
set of states as the exact original Hamiltonian ||. From among the eigenstates of 
the Hamiltonian (3), it is necessary to choose only those that correspond to the 
same center-of-mass energy. This can be achieved first by working in a complete 
Nhtt basis space, and then by projecting the center-of-mass eigenstates with en¬ 
ergies greater than | hfl (representing spurious center-of-mass motion) upwards 
in the energy spectrum In our case, we do this by adding (3PH^} m P to 
and subtracting fi^frPlP from equation (3) above. One unit of P,p m has already 
been incorporated, as mentioned above, and it is also subtracted at this stage. 
The resulting shell-model Hamiltonian takes the form 


= P H 

i<j 
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where /3 is a sufficiently large positive parameter. 

When applied in a suitably chosen finite HO basis space, this procedure 
removes the spurious center-of-mass motion, and has no effect on the intrinsic 
spectrum of states ||. Our choice of the P-space is fixed by NM1 that signifies 
a many-body basis space of all HO Slater determinants having total HO quanta 
less than or equal to the cutoff given by N. We refer to this as a complete Nhfl 
basis space. It is complete in the sense it allows the precise removal of spurious 
center-of-mass motion effects. 

In principle, the effective interaction introduced in Eqs. (3) and (4) above 
should reproduce exactly the full-space results in the finite basis space for some 
subset of physical states. Furthermore, an A —body effective interaction is re¬ 
quired for an A —nucleon system. In practice, however, the A —body effective 
interaction cannot be calculated exactly, and it is approximated by a two-body 
effective interaction determined for a two-nucleon sub-system or two-nucleon 
cluster. 

In this work, we follow the procedure described in Refs. H in order to 
construct the two-body effective interaction. The procedure employs the Lee- 
Suzuki 111 similarity transformation method, which yields an interaction in the 
form 


P 2 V eff P 2 = P 2 VP 2 + P 2 VQ 2 uP 2 , 

with oj the transformation operator satisfying uj = Q 2 u)P 2l and P 2 , Q 2 = 1 — 
P 2 operators that project on the two-nucleon finite basis and complementary 
spaces, respectively. Note that we distinguish between the two-nucleon system 
projection operators P 2 , Q 2 and the A—nucleon system projection operators 
P, Q. Note also that the size of the two-nucleon basis is dictated by the choice 
of Nhfl for the A-body problem. 

Our calculations start with exact solutions of the two-body Hamiltonian 


Hh‘ = //, 


02 


Pi 


Vo 


■P2 


(5) 


H—+ r 
2m 2 v ; 


Vi 2 (fi - r 2 ) - - r 2 ) 2 


which follows from reducing the shell-model Hamiltonian (2) to the two-nucleon 
case. The exact solutions are obtained by assuming a OhQ motion of the center- 
of-mass of nucleons 1 and 2. Let | ap) denote the two-nucleon HO states which 
form the basis space signified by P 2l and let \o.q) denote those that form the 
Q 2 — space. Dropping the subscript “2” for convenience, one can then express 
the Q —space components of an eigenvector | k) of the Hamiltonian (5) as a 
combination of the P—space components with the help of the operator u: 


{a Q \k} = ^2(a Q \uj\ap)(ap\k) 

a P 


( 6 ) 
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If the dimension of the basis space is dp , we may then choose a set K oi dp 
eigenvectors for which relation (6) will be exactly satisfied |I2|. Typically, these 
dp states will be the lowest states obtained in a given channel. Under the 
condition that the dp x dp matrix (ap\k) for \k) e K is invertible, the operator 
uj can be determined from Eq. (7), and the effective Hamiltonian can then be 
constructed as follows: 


{lp\H 2 eff\ap) = ^2 

k e K 


{'yp\k)Ek(k\ap) 


yZ('yp\k)Ek(k\aQ}(aQ\uj\ap) 

OLQ 


(7) 


This Hamiltonian, when diagonalized in a finite basis space, reproduces exactly 
the set K oi dp eigenvalues Ek of the defined 2-body problem. Note that the 
effective Hamiltonian (7) is not Hermitian, and that it can be transformed into a 
Hermitian form H 2e f / by applying a similarity transformation determined from 
the metric operator P 2 (l + uj^lo)P 2 ju|: 

Hieff = [P 2 (l +^u)P 2 } 1/2 H 2eff [P 2 [l + ^co)P 2 }~ 1/2 

The two-body effective interaction actually used in the present calculations is 
determined from this two-nucleon effective Hamiltonian as V 2e f / = H 2e f f — Hq 2 ■ 
The resulting two-body effective interaction V 2 e // depends on A, on the HO 
frequency H, and on N max , the maximum many-body HO excitation energy 
(above the lowest configuration) defining the P— space. Furthermore, when used 
in the shell-model Hamiltonian (4), it results in the factorization of our many- 
body wave function into a product of a center-of-mass | hfl component times an 
intrinsic component, which allows exact correction of any observable for spurious 
center-of-mass effects, thus preserving translational invariance. This feature 
distinguishes our approach from most phenomenological shell-model studies that 
involve multiple HO shells. 

So far, the most important approximation used in our approach is the neglect 
of contributions coming from higher than two-body clusters to our effective 
Hamiltonian. Although our method can be readily generalized to include the 
effects of three-body clusters, leading to the derivation of a three-body effective 
interaction, so far this has only been done in A = 3 and A = 4 systems 00. 
and it is not yet clear how difficult it will be to extend it to the p— shell nuclei 
with A > 5 and with the same range of NhQ values accessible at the 2-body 
cluster level. 

In order to approximately account for the many-body effects neglected when 
using only a two-body effective interaction, we introduce a simple residual inter¬ 
action that depends on all particles. More specifically, our residual many-body 
interaction contains a delta function times the sum of the oscillator quanta for 
the A — 2 spectators denoted by N sps and defined by 
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N sps = N sum -N a = N' sum - N-y (8) 

where N sum and N' sum are the total oscillator quanta in the initial and final 
many-body states, respectively. The quantities N a and TV 7 are the total oscil¬ 
lator quanta in the initial and final two-nucleon states |a) and I 7 ), respectively, 
interacting through the delta function. 

Why do we introduce such a dependence of the many-body forces? We are 
motivated by the experience with the “multi-G” method of addressing effective 
many-body forces Jf|, which is based on the role of the spectators in a many- 
body matrix element. As the spectators are excited, the effective P 2 space for 
the interacting pair decreases and this increases the attraction of the G-matrix. 
The spectator dependence in the multi-G approach to the effective many-body 
forces predominantly lowers the excited multi-?ifl states. 

In this connection, one might imagine, as we do, the spectator dependence 
to be of a form where the strength is weakest when the spectators have Ohfl 
excitation (above their minimum value) and strongest when they have the max¬ 
imum excitation. These two situations correspond, respectively, to when the 
2-particles roam over the entire NhSl and over only the lowest hfl range. Thus, 
we expect this 7V sps -dependent term to be attractive and account for the leading 
effects observed in the multi-G investigations. 

We also allow for a spectator-independent term in the many-body force 
and we will take it to be a delta function for simplicity. We expect this term 
will approximate the many-body effective force components that are largely 
independent of the spectator configurations. We expect this contribution to be 
repulsive in order to correct for the overbinding found in smaller model spaces 
when H e f / is approximated with the 2-body cluster form ||. We will see below 
that the fits to nuclear spectra verify these expectations. 

Our total effective many-body interaction is then expressed in the form of 
a zero-range two-body interaction, and two parameters, one of which multiplies 
the many-body variable N sps introduced above: 

Vres{N sps ) = (a + bN sps ) x GS{f{ - 7=2) ( 9 ) 

where G is a fixed strength parameter taken to be 1 MeV fm 3 in the present 
work. The parameters a and b are basis-space dependent dimensionless param¬ 
eters to be determined by fits to selected many-body observables. 

We will then use this schematic residual interaction to study the role played 
by the many-body forces that were neglected in our derivation of the effective 
two-body interaction. We shall derive the latter from the realistic Argonne 
VS' NN potential ||Tg|, and perform calculations of some basic observables in 
the A = 4 system with and without the addition of the residual many-body 
interaction V res (N sps ). and in various basis spaces. We solve for the properties 
of A = 4 nuclei using the m— scheme Many-Fermion Dynamics code [ [iT] that 
was developed for No-Core Shell-Model (. NCSM ) calculations. 
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3 Application to A = 4 Nuclei 

In this section we apply the methods outlined in Sec. II to derive the effective 
two-body interaction, based on the Argonne V8' NN potential, which we then 
use in the remainder of our work. However, it is now well-known that two-body 
interactions alone are not sufficient, and in order to better account for such 
experimental facts as the binding energy of nuclei, one is led to introduce three- 
body and higher multi-body forces. In the following discussion, we first consider 
calculations done with only the Argonne V8 1 two-body interaction in order to 
see how such two-body calculations compare with experiment. We then ask the 
following question: what couplings (a, b) for the residual many-body interaction 
are needed in order to achieve as close an agreement with experiment as possible? 

We use a complete NhQ basis space with, e.g., N = 8 for the positive- 
parity states of 4 He. This means that nine major HO shells are employed. 
The two-nucleon basis space is defined in our calculations by N max , e.g., for 
an 8 htt calculation for 4 He we have N max =8 . The restriction of the HO shell 
occupation is given by (A/ + N 2 ) < N max where the single-particle HO quanta 
are given by W = 2m + The same conditions hold for the relative 2 n + l, 
center-of-mass 2N + L , and 2 n + l + 2N + L = Ni + N 2 quantum numbers, 
respectively. 

In the case of the A = 4 system, we seek to reproduce the following oberv- 
ables as accurately as possible: the binding energy of the ground state ( GS ) as 
well as the energies of the first three excited states. We are also interested in 
reproducing the root-mean-square radius R. 

3.1 The Experimental Situation 

In Table 1 we present the experimental energies of the low-lying states that are 
relevant to our present study. The binding energies of the ground states are 
given first, followed by the excitation energies of the first three excited states 
(all energies are in MeV) 0 0 - Furthermore, the root-mean-square point- 
charge radius of the ground state of 4 He is known experimentally to be 1.46 
/to, which we also take to be the point mass radius. 

3.1.1 Calculations Using Only Two-Body Effective Interactions 

Before we embark upon the investigation of the many-body force, let us first 
examine the theoretical predictions of the 4 He spectrum using only the Argonne 
V8' plus the Coulomb interaction. In Table 2, we present results of calculations 
of the first two positive-parity and the first two negative-parity states, using the 
interaction in increasingly larger basis spaces (indicated by N rnax . as defined 
above). The results depend somewhat on the choice of hQ, the HO energy 
parameter, but we shall adopt the value Ti fl=19 MeV , which is suggested by 
the phenomenological expression Q 
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In the large N max limit our results become independent of the choice of TiVL j2lj . 

As shown in Table 2, it is clear that the Argonne V8' interaction overbinds 
the ground state for N max < 8 and underbinds it for N max = 8 and 10, with the 
overbinding decreasing (and the underbinding increasing) as N max increases. 
On the other hand, the excitation energy of the 0^" state steadily decreases 
with N max . If a two-body force is alone sufficient, and if the Argonne 148' plus 
Coulomb interaction represents a sufficiently accurate form of this two-body 
force, we would expect our calculations to converge to the experimental val¬ 
ues. By examining Table 2, it becomes clear that while the first positive-parity 
excited state (0^) and the negative-parity states (0 _ and 2“) are becoming 
more bound with increasing Nmax, the ground state itself ( GS ) is becoming less 
bound. 

Furthermore, while the discrepancy between the calculated and the exper¬ 
imental ground-state energy is reduced with increasing N max (from 12% with 
N ma x = 2 to less than 4% with N max = 10), the improvement is less favorable 
for the energies of the excited states ( e.g ., from 51% with N max = 2 to 23% with 
Nmax = 10 for the 0^" state). The percentages used here represent the percent 
error (in absolute value) of the calculated energy relative to the experimental 
value. For the ground state, this is based on the total binding energy, whereas 
for the excited states it is based on the excitation energy. 

We are, thus, led to conclude that the two-body force used is somewhat 
deficient and we need either a better two-body force, or three-body and/or 
higher many-body forces, in general. This is exactly the situation that we wish to 
investigate with our zero-range, spectator-dependent, residual delta interaction. 


3.1.2 The effect of adding a residual many-body force 

Next, we consider the effect of adding our residual delta interaction to the de¬ 
rived effective 2-body interaction in order to obtain a 2 + many-body potential: 

V = V A rg 0nne V8' + VCoul + [a + bN sps ] X GS(f{ — f^). 

By adding the spectator-dependent phenomenological interaction introduced in 
Eq. (9), we are simply modelling the contribution of all residual many-body 
forces, omitted by our derived two-body effective interaction. We then vary the 
coefficients a and b independently for each basis space, in order to achieve good 
agreement with the experimental ground-state binding energy and excitation 
energy of the first excited 0 + state of 4 He. The results are shown in Table 3 
and illustrated in Fig. 1. We note that a is positive, while b is negative, and that 
both are decreasing functions of N m ax, with a decreasing somewhat faster than 
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b. As the ground state was found to be overbound, whereas the first excited 
state was found to be underbound by the pure two-body interaction (see Table 
2 ), it should not be surprising that a turns out to be positive, while b is negative. 
It is also reassuring to see both a and b decrease steadily (in absolute value) 
with N max , indicating that the need for the residual many-body interaction V res 
decreases and that the adequacy of the two-body interaction increases as the 
basis space becomes larger. More precisely, the result that a —> 0 as N max — > oo, 
is not surprising, as it indicates that there is little or no correction needed at 
the pure two-body level to our Hamiltonian. On the other hand, the 6 -term, 
which carries all the information about the ’’true” many-body interaction, is not 
converging to zero. We, therefore, conclude that it is possible to accommodate 
the contribution of the missing many-body interaction and to represent it in the 
simple form of bN sps x GS(f{ — r 2 ) for these observables. 

We are now in a position to examine the predictive power of our extended 
effective Hamiltonian. We begin by asking how well such a residual many-body 
interaction (determined in the case of 4 He) can help reconcile the calculated 
spectra of 4 H and 4 Li with the experimental situation? In Tables 4 and 5 we 
show the results obtained with and without the residual many-body interaction 
and compare them to experiment. We are now addressing (T = 1) negative- 
parity states, so we perform our calculations in an N max = 9 basis space, and 
we use the values a = 0.63 and b = —2.43, which were determined in the earlier 
fit to 4 He (see Table 3). 

It is clear from the results shown in Tables 4 and 5 that a much better 
agreement with experiment is achieved for both 4 H and 4 Li when our residual 
many-body interaction is added to the pure two-body interaction. Aside from 
the improved agreement for the ground states and first excited states, we ob¬ 
serve that it is only by including the residual many-body interaction that the 
calculated 0~ and 1^ states in 4 Li become bound, in accord with experiment. 

Based on the overall results for the A = 4 isotopes, presented above, we 
can identify a procedure for improving the many-body force. In particular, the 
4 He fit can be interpreted as determining the spin-isospin averaged values of 
the couplings for our many-body force. By allowing spin-isospin dependence in 
the couplings, we could fine tune the coupling strengths, while preserving their 
average values, so as to obtain an improved description of the other isotopes. 
Thus, we can introduce an isospin dependence in our residual many-body force 
as follows: 


V res = [(a + a'fi ■ t 2 ) + (6 + b'n ■ f^)IV S p S ] x G5(f[ - rj). 

or, equivalently, by assigning different values for a and b in the T = 0 and T = 1 
channels, denoted by a 0 , a\ and bo, b±. We do this while keeping the isospin 
average of ao and ai to be equal to a, and likewise for bo and b 1 . Applying 
a minimum root-mean-square convergence criterion on the calculated energies 
of the first two states in each of the A = 4 nuclei considered relative to their 
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experimental values, we find that it is possible to improve the calculations when 
using an isospin-dependent version of the residual many-body interaction (Table 
6 ). 

It would be interesting to study whether our model of the residual many- 
body interaction can be generalized to the study of A > 4 systems. For example, 
the test of the spin-isospin dependent many-body force can then be carried out 
in the A = 5 systems. This effort is underway and will be reported in another 
paper. 

To summarize, we have presented a simple model for treating multi-body 
interactions in nuclear-structure calculations in terms of a two-parameter, basis- 
space-dependent delta interaction. We determine the two-parameters as a func¬ 
tion of basis-space size by fits to the lowest two 0 + states of 4 He. The parameter 
values decrease with basis-space size, as anticipated. The obtained multi-body 
interaction also yields improved agreement with experiment for negative parity 
states in 4 He and for the T = 1 states of 4 H and 4 Li, which were not included 
in the empirical fits. The advantages of our model for the multi-body interac¬ 
tion in nuclei are two-fold: (1) the model provides a simple method for including 
multi-body interactions into nuclear-structure calculations, thereby yielding im¬ 
proved agreement with experiment, and (2) the model will give us insight into 
the nature of such multi-body interactions, so that we will be able to model 
them more realistically in the future. 
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Figure Captions 

Fig. 1 : The strengths a (solid curve) and b (dashed curve) of the residual 
delta interaction needed to achieve agreement with experiment in the case of 
4 He and using hfl = 19 MeV. 
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Table 1: Experimental binding energies of the ground states ( GS ) and excitation 


energies of the first three excited states (in MeV) in the A = 4 system. 


4 * 
II ^ 
o 

3 

(T 

Li 

= i) 

(T = l) 

0+ (GS) -28.295674 5 

2~ (GS) 

-4.618058 

2b (GS) -5.575154 

0 ^ 20.210 

ir 

0.320 

lb 0.310 

0 b 21.010 

0 b 

2.080 

0 b 2.080 

2(" 21.840 

1-2 

2.850 

lb 2.830 


Table 2: Ground-State energy and excitation energies (in MeV) of the first two 
0 + states and the first two negative-parity states of 4 He using only the Argonne 
V8 1 and the Coulomb potentials ( i.ea — 0 and b = 0), as calculated in shell- 
rnodel spaces characterized by N max htt excitations and with Ml = 19 MeV. 


(a) 

Positive Parity 

(b) Negative Parity 

N 

1 v max 

E(GS) 

E x (0t) 

L*max 

E x ( 0") 

E x ( 2“) 

2 

-31.853 

30.577 

3 

26.560 

28.692 

4 

-30.758 

29.825 

5 

25.694 

27.600 

6 

-29.328 

27.019 

7 

23.884 

25.644 

8 

-28.149 

25.213 

9 

22.607 

24.211 

10 

-27.247 

24.902 
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Table 3: Values of the parameters a and b in the interaction V^J gonneV8 , + 
Vcoui + [a + b * N sps ] x G5(ri — rj) and corresponding energies as calculated in 
various basis spaces characterized by N max hfl excitations, with hQ = 19 MeV 
for 4 He. _ 


a 

b 

(a) Positive Parity 

(b) Negative Parity 


N 

iy max 

E(GS) 

Rgs 

E x { 0J) 

N 

iy max 

E x ( 0") 

E x { 2“) 

2.22 

-6.86 

2 

-28.305 

1.390 

20.223 

3 

23.811 

26.017 

2.09 

-5.80 

4 

-28.291 

1.341 

20.221 

5 

23.638 

25.353 

1.15 

-3.12 

6 

-28.296 

1.354 

20.212 

7 

22.446 

23.977 

0.63 

-2.43 

8 

-28.301 

1.383 

20.238 

9 

21.541 

22.903 

0.09 

-1.78 

10 

-28.295 

1.409 

20.209 





Table 4: Experimental and calculated energies (in MeV) of the first few states in 
4 77, with and without the residual many-body interaction (a = 0.63, b = —2.43) 
in an N max = 9 basis space and using Ml = 19 MeV. 


State 

Experimental 

Energy 

Energy calculated without 
residual many-body interaction 

Energy calculated with 
residual many-body interaction 

2r {GS) 

-5.575 

-2.994 

-4.922 

ir 

-5.265 

-2.452 

-4.792 

or 

-3.495 

-1.403 

-3.861 

1-2 

-2.745 

-0.935 

-3.864 


Table 5: Experimental and calculated energies (in MeV ) of the first few states in 
4 Li, with and without the residual many-body interaction (a = 0.63, b = —2.43) 
in an N max = 9 basis space, and using HQ = 19 MeV. _ 


State 

Experimental 

Energy 

Energy calculated without 
residual many-body interaction 

Energy calculated with 
residual many-body interaction 

2r (GS) 

-4.618 

-1.409 

-3.433 

ir 

-4.298 

-0.933 

-3.343 

or 

-2.538 

0.097 

-2.422 

1-2 

-1.768 

0.532 

-2.444 
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Table 6: Experimental and calculated energies (in MeV) of the first two states 
in the A = 4 system, with the isospin-independent and the isospin-dependent 
residual many-body interaction in N max = 8, 9 basis spaces, and using htt = 


19 MeV. 


Nucleus (State) 

Experimental 

Energy 

Using ao = ai = 0.63 
and bo = bi = —2.43 

Using oo = 0.52, ai = 0.67 
and bo = —2.88, bi = —2.28 

4 He(0+ ( GS )) 

-28.297 

-28.301 

-28.627 

A He{ 0+) 

-8.086 

-8.063 

-8.697 

A H{ 2r (GS)) 

-5.575 

-4.922 

-5.255 

4 H( ir) 

-5.265 

-4.792 

-5.157 

4 Li(2i (GS)) 

-4.618 

-3.433 

-3.775 

ir) 

-4.298 

-3.343 

-3.713 
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